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ABSTRACT 

As part of a large Spitzer/MIPS GTO program, we have searched for infrared 
excesses due to debris disks toward 26 FGK field stars known from radial velocity 
(RV) studies to have one or more planets. While none of these stars show excesses 
at 24 /im, we have detected 70 /im excesses around 6 stars at the 3-a confidence 
level. The excesses are produced by cool material (< 100 K) located beyond 10 
AU, well outside the "habitable zones" of these systems and consistent with the 
presence of Kuiper Belt analogues with ~100 times more emitting surface area 
than in our own planetary system. These planet-bearing stars are, by selection 
for RV studies, typically older than 1 Gyr, and the stars identified here with 
excesses have a median age of 4 Gyr. We find a preliminary correlation of both 
the frequency and the magnitude of dust emission with the presence of known 
planets. These are the first stars outside the solar system identified as having 
both well-confirmed planetary systems and well-confirmed IR excesses. 

Subject headings: infrared:stars; Kuiper Belt; planetary systems:formation; plan- 
etary systems:protoplanetary disks 



1. Introduction 

First discovered by the Infrared Astronomical Satellite (IRAS), the "Vega phenomenon" 
is characterized by a significant deviation from the expected Rayleigh- Jeans emission of 
quiescent main-sequence stars (Aumann & Good 1990; Gillett 1986). The excess IR emission 
from other stars has been attributed to dust orbiting at distances anywhere from less than 1 



- 2 - 



to more than 100 AU, analogous to the zodiacal cloud and the Kuiper belt in our own solar 
system. Due to effects such as radiation pressure and Poynting-Robertson drag, the lifetime 
of the dust is generally much shorter than the age of the system; any dust observed must 
have been recently produced. In the solar system, dust is generated by collisions between 
larger bodies in the asteroid and Kuiper belts, as well as from outgassing comets. Extrasolar 
systems with IR excess presumably have their own supply of large, solid planetesimals and, 
perhaps, large planets like those in our system. 

IRAS found that approximately 15% of main-sequence stars (PIcts & Vynckier 1999), 
predominantly A and F stars, show excesses attributable to dust in orbit around the stars 
at distances between a few and a few hundred AU and heated by absorbed starlight to 
temperatures of 30-300 K (Backman & Paresce 1993). Fractional luminosities, Ldust/-^'*; are 
observed as large as 10~^ for stars with bright excesses such as /3 Pic and as small as 10"^ 
for stars at the detection threshold. These values can be compared to our solar system: 10~^ 
to 10~^ for the hot material in the zodiacal cloud (1-3 AU; Dermott et al. 2002) and 10~^ 
to 10~^ predicted for the cooler Kuiper Belt region beyond 10 AU (Stern 1996), which is 
constrained to be < 10^^ by IRAS (Aumann & Good 1990; Tcplitz et al. 1999). Observations 
by the Infrared Space Observatory (ISO; Habing et al. 2001; Spangler et al. 2001) suggest a 
decline in the fraction of stars with excess IR emission with time, but with the possibility 
of finding modest (Ljust/-^* > 10"^) excesses among older stars (Decin et al. 2000, 2003). 
Spitzer and IRAS results toward more than 200 A stars confirm a general decline in the 
average amount of emission with little 24 /im emission seen for sources older than 100 Myr 
but with large variations within each age group studied (Rieke et al. 2004). The variations 
arc probably due in part to sporadic replenishment of dust clouds by individual collisions 
between large, solid bodies, but also are likely to reflect a range in mass and extent for the 
initial planetesimal disk. 

Images in the optical and near-IR (Heap et al. 2000; Schneider et al. 1999), submillime- 
ter (Holland et al. 1998; Greaves et al. 1998; Holland et al. 2003) and from Spitzer itself 
(Stapelfeldt et al. 2004) have revealed that in some cases the dust exists in flattened disks 
with warps, gaps, and blob-like structures that have been attributed to the gravitational 
influence of planets (Wyatt & Dent 2002; Dermott et al. 2002). Still, while the frequency 
and magnitude of IR excess are strongly correlated with stellar youth, no such correlation 
has yet been found between IR excess and the presence of planets. Even before the first 
confirmed discovery of an extrasolar planet orbiting a Sun-like star (51 Peg b; Mayor & 
Queloz 1995; Marcy & Butler 1995), attempts had been made to find links between planets 
and dusty debris around other stars (e.g. e Eri; Backman & Paresce 1993). While early 
radial-velocity searches for planets found a planet-like trend for e Eridani (Gampbell et al. 
1988), the candidate planet remains controversial even today due to the extremely high level 
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of chromospheric activity on this young star. The 7-year periodicity found in the star's radial 
velocity (Campbell et al. 1988; Gumming et al. 1999; Hatzes et al. 2000) may be induced by a 
planet orbiting at several AU (Hatzes et al. 2000) or may instead be a reflection of the star's 
poorly understood magnetic cycle. Conversely, 55 Cnc has a well-known planetary system, 
but claims of a dust disk based on visible imaging (Trilling et al. 2000) and ISO photometry 
(Dominik et al. 1998) have been ruled out by HST coronography (Schneider et al. 2001), sub- 
millimeter mapping (Jayawardhana et al. 2002), and Spitzer/MlPS measurements presented 
here. In spite of the many studies specifically designed to find a connection between dust 
and planets, no systems with both well-confirmed planets and well-confirmed IR excess have 
yet been observed. In fact, the data to date appear to indicate an anti-correlation between 
dust and planets (Greaves et al. 2004). 

The longstanding assumption that stars besides the Sun might harbor both planets 
and extended dust emission has still not been verified. The aim of this paper is to use the 
Multiband Imaging Photometer on the Spitzer Space Telescope (MIPS; Rieke et al. 2004) 
to directly address this outstanding issue by searching for excess IR emission from known 
planet-bearing stars. MIPS provides unprecedented sensitivity at mid-IR wavelengths and 
as such is an ideal instrument for searching for the dust emission from solar-type stars. 
While the focus of this paper is on observations of stars which are known to have planets, 
for purposes of noise calibration and to greatly improve our statistical analysis, the results 
for planet-bearing stars are considered within the context of a larger program - the Volume 
Limited Survey (VLS) . VLS is a Spitzer GTO program designed to search for excesses around 
a broad sample of 155 nearby, F5-K5 main-sequence field stars, sampling wavelengths from 
8-40 //m with IRS and 24 & 70 //m with MIPS. This survey consists of two components 
- stars with planets and those without. The stars without known planets represent a true 
volume-limited sample selected for expected 70 /im brightness and low levels of infrared 
cirrus. The sample includes G and late F stars out to 25 pc and early K stars out to 15 pc. 
This component of VLS consists of 114 stars (without planets), 58 of which have already 
been observed. These stars will be discussed in more detail in a separate paper (Bryden et 
al. 2005); here we focus on the second component of the VLS program - planet-bearing stars. 
The planet-bearing sample was selected from a group of stars of similar spectral types known 
to have RV planets. Of 41 planned targets, 26 have been observed to date. These stars are 
not necessarily located in regions of low infrared cirrus and are typically farther away (~15- 
50 pc) than the non-planet VLS sample. It is worth noting that the disk around e Eri, a KO 
star located 3.3 pc away, would have been undetectable by IRAS if that star were located 
even 10 pc away; its 60 /xm flux density would then have dropped below the completeness 
hmit of the Faint Source Gatalog (Moshir et al. 1990). The sensitivity of Spitzer/MIPS 
allows the search for disks not just around the nearest FGK stars, but also around those 
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located at distances of 10-50 pc around which planets have been identified. 

We present our MIPS observation in §2, along with a full error analysis to determine 
whether any measured excesses are statistically significant. For the systems with IR excess, 
§3 attempts to find a correlation with system parameters such as stellar metallicity/age 
or planet mass/semi- major axis. Before the concluding remarks (§5), we describe how our 
MIPS observations constrain the dust properties in each system (§4). 

2. Observations 

During its first year of operation Spitzer observed from the VLS sample 26 stars with 
RV planets (Table 1) and an additional 58 stars without any known planets. All stars were 
observed at 70 //m. One planet-bearing star, HD128311, was not observed at 24 //m due to 
scheduling constraints, but its spectrum at this wavelength was measured with IRS Long- 
Low and was found to be consistent with no excess shortward of 30 /xm. The remainder of 
the planet-bearing stars were observed at 24 jj,m and detected with high signal-to-noise ratio 
(SNR). 

2.1. Data Reduction 

Our data reduction is based on the DAT software developed by the MIPS instrument 
team (Gordon et al. 2004a). Because the data pipeline is still under development and the 
adopted calibration numbers are likely to change with further improvements in the overall 
MIPS data analysis, we list our chosen set of data reduction parameters in detail. 

For each of the 24 fim and 70 bands, the instrument calibration was derived from 
~75 observations of ~35 different stars. To minimize systematic uncertainties, a mix of 
stellar types was used: most of the calibrators are K giants, several are A stars, and a few 
are solar analogs. The 24 /xm photosphere predictions were derived from stellar templates 
(e.g. Cohen et al. 1999) for K giants, from Kurucz models (e.g. Kurucz 2003) for A stars, 
and from the solar spectrum (e.g. Tobiska et al. 2000) for the solar analogs. The 70 //m 
photospheric predictions were obtained similarly, except that the analytic spectral functions 
of Engelke (1992) were used to extrapolate the K giant templates out to 70 /xm. 

At 24 yum, we carried out aperture photometry on reduced images using an aperture 6 
camera pixels (1 pixel=2.5" at 24 /im) in radius, a background annulus from 12 to 17 pixels, 
and an aperture correction of 1.15. The fiux level is calibrated at 1.047 //Jy/arcsec^/(DN/s), 
with a default color correction of 1.00 appropriate for sources warmer than 4000 K at a 
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weighted-average wavelength of 23.68 /xm. Final images were mosaiced from individual 
frames with quarter-pixel subsampling. 

At 70 iim we used images processed beyond the standard DAT software in order to cor- 
rect for time-dependent transients, corrections which can significantly improve the sensitivity 
of the measurements (Gordon et al. 2004b). At 24 //m, a large enough aperture was chosen 

so as to integrate all of the flux out to roughly the first Airy ring. Because the accuracy of 
the 70 fim data is limited by background noise, rather than instrumental effects, such a large 
aperture contains an undesirable amount of sky fluctuations. In order to maximize signal-to- 
noise, a smaller aperture is more appropriate at 70 fim. - just 1.5 pixels in radius. With a 4 to 
8 pixel radius sky annulus, this aperture size requires a relatively high aperture correction of 
1.79. The flux level is calibrated at 15,800 /iJy/arcsec^/MIPS_70_unit, again with a default 
color correction of 1.00 (MIPS_70_units are based on the ratio of the measured signal to the 
stimulator flash signal). Images were mosaiced from individual frames with half- pixel sub- 
sampling. For both the 24 //m and 70 //m data, neighboring point sources were subtracted 
from the images before measuring the sky brightness. Also for both wavelengths, we fix 
each target's position, rather than centroid onto the peak of the local flux distribution. If 
the target position is not flxed, the centroid position tends to float toward any local noise 
fluctuations, resulting in an inappropriately high flux measurement, particularly for noisy 
fields. The fact that no centering of the aperture is required is a refiection of the excellent 
telescope pointing intrinsic to Spitzer (<1"; Werner et al. 2004). 

In order to determine whether any of our target stars have an IR excess, we compare 
the measured photometry to predicted photospheric levels. For convenience in extrapolating 
the photospheric emission we have fitted Kurucz models (CasteUi 2003; Kurucz 2003) to a 
combination of short wavelength data taken primarily from the Tycho/Hipparcos (visible), 
Johnson (bright near-IR) or 2MASS (faint near-IR), and IRAS (12 //m) catalogs. After 
excluding one outlying point for a star with an excess (HD69830, which will be discussed 
elsewhere; Beichman et al. 2005), the 83 flux measurements at 24 //m have an average 
-^MiPS /-^photosphere of 0.99 which, although gratifyingly close to unity, is not surprising given 
that the present Spitzer calibration is based on similar stellar models. More importantly for 
our purposes, the dispersion of -Fmips /-^photosphere is 0.07 for this sample (Fig. 1). 

At 70 iim, 73 out of 84 stars (both planet-bearing and without planets) are detected 
with signal-to- noise ratio > 3. This is in contrast with previous IR surveys of AFGK stars 
with ISO, in which only half of the stars without excess were detected (Habing et al. 2001). 
While the sensitivity of these Spitzer observations is roughly a factor of 10 better than 
previous data, Spitzer's sensitivity is limited by extragalactic source confusion and cirrus 
(see below) which will make it difficult to look for excesses around stars much fainter than 
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those discussed here. 

The distribution of 70 pLm flux densities relative to the expected photospheric values is 
shown in Fig. 2. Unhke the tight distribution of flux ratios at 24 //m, several stars have 70 
/im flux much higher than expected from the stellar photosphere alone. In one case, the 70 
//m flux is high by a factor of 15. Many of the stars (12) will be identifled in the following as 
having statistically significant IR excess. Ignoring these stars with excesses and those with 
SNR < 3, the average ratio of MIPS flux to predicted photosphere is |-Fmips /-^photosphere I = 
1.03 ± 0.04, consistent with the present calibration. More important than the absolute 
calibration of -Fmips /-^photosphere is its dispersion, which indicates the overall accuracy of the 
flux measurements and gives an immediate sense of the minimum level of excess necessary 
for detection. The dispersion in the 70 iim sample is 28% (excluding the stars with excesses), 
considerably higher than that in the 24 //m data (7%). The following section (§2.2) discusses 
in more detail the noise levels within the 70 /im data. 

Fig. 3 shows the spectral energy distributions for two of the planet-bearing stars, 
HD82943 and HDl 17176, that have excess 70 //m emission. Published photometric fluxes 
for these stars from visible to infrared are well flt by Kurucz stellar atmospheres (dotted 
fines; Castelfi 2003; Kurucz 2003). The Spitzer/MlPS 24 pim fluxes are also weU flt by the 
model atmospheres (within 1% for HD82943, 5% low for HD117176). Tfie 70 /xm emission, 
meanwhile, is well above that expected from the stellar photosphere alone; an additional 
component of emission due to dust must be added to the Kurucz model in order to fit the 
70 fim data. Because there is only a single 70 fim. measurement of IR excess, however, the 
SED can be fit by a range of dust temperatures and luminosities. Constraints on the dust 
in these systems are discussed in §4. 



2.2. Detection of Excess at 70 /xm 

A potential correlation of planets with IR excess is already apparent in Fig. 2, in which 
4 of the 5 largest 70 fim fiux ratios are for stars which harbor planets, even though planet- 
bearing stars make up less than one-third of the sample. Still, without a proper analysis of 
the noise levels in each field it is impossible to assess whether the IR excesses are statistically 
significant. 

In general, the sources of noise can be divided into pixel-to-pixel noise (which we refer 
to as "sky" noise) and errors for the field as a whole ( "systematic" noise) . The pixel-to-pixel 
noise is a combination of detector/photon noise along with real sky background fiuctuations. 
For the 24 //m measurements, photon noise is negligible. Even with the minimum allowed 
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integration time (1 cycle of 3 sec exposures = 48 sec), the sensitivity of MIPS is overwhelming; 
even our dimmest source could theoretically be detected in just a few milliseconds. Further, 
the background noise is low: galactic cirrus is weak at this wavelength; the zodiacal emission 
is relatively smooth across the field of view; and the confusion limit for distant extragalactic 
sources is just 0.056 mJy (Dole et al. 2004b). 

Instead, for the 24 //m measurements, systematic errors dominate. The instrumental 
contribution to these errors is thought to be very low: 24 yum observations of calibrator stars 
are stable with 1% rms deviations over several months of observations (Rieke et al. 2004). 
In addition to any residual instrumental problems, the dispersion in Fmips /-^photosphere can 
include errors in the photosphere extrapolation and the effects of source variability. The 
fitting of the photosphere can be as precise as 2% when the best 2MASS photometry is 
available (Skrutskie et al. 2000), but for stars brighter than < 6 mag, 2MASS data are less 
accurate and lower precision near-IR data and/or shorter wavelength observations must be 
used. Extrapolation from visible data places considerably more weight on the photospheric 
models and may result in greater uncertainty in the predicted photospheric levels at 24 /im 
and 70 /im. 

While the relative contributions of instrumental and fitting errors to the 24 /im disper- 
sion are unknown, it is clear that the photometry is currently limited to 7% accuracy due 
to their combined effects. While this may improve with further instrument analysis and 
better photospheric models, 3-(T detections of excess at 24 fim now require measured fluxes 
at least 20% above the stellar photosphere. Only one of our sources (HD69830) meets this 
requirement, with a strongly significant 45% excess. The relatively hot dust emission in this 
system, which is not yet known to have any planets, will be discussed in a separate paper 
(Beichman et al. 2005). 

At 70 fim we planned the number of observing cycles (each with 126 sec integration time) 
so that the expected detector /photon noise would be ~0.1 of the stellar fiux, or equivalently, 
to achieve a signal-to-noise (SNR) ~ 10 on the stellar photosphere. Integration times range 
from a few minutes up to 1 hour. Within each reduced image, we estimate the sky noise 
by directly measuring the standard deviation in the background flux when convolved with 
our chosen aperture size. Based on this measured sky noise, wc flnd a median SNR of ~6, 
excluding the sources identifled as having an excess. This noise estimate, however, includes 
contributions beyond just detector/photon noise. Unlike at 24 /im, the fluctuations in the 
background at 70 fim can be signiflcant. This background, a combination of galactic cirrus 
and extragalactic confusion, creates a noise floor for each fleld which cannot be improved 
with increased integration time. To try to minimize this problem, the target stars in the 
Volume Limited Sample came from areas of low galactic cirrus, as estimated from the IRAS 
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Sky Survey Atlas (IPAC 1994). The confusion limit for extragalactic background sources, 
however, is unavoidable and sets a strict lower limit for the sky noise at 70 fira. Dole et al. 
(2004b) find a 5-a confusion limit of 3.2 mJy based on source counts from Spitzer (Dole 
et al. 2004a) extrapolated to fainter fiuxes with the model of Lagache et al. (2004). In 
our sample, the lowest (l-cr) noise levels observed toward stars located in clean portions of 
the sky are ~2 mJy, somewhat worse than Dole et al.'s best-case limit. This difference is 
attributable to the larger effective beam size used for our photometry/noise calculations. 
On top of the confusion limit, a few sources have higher noise values due to galactic cirrus 
and/or modestly worse than typical detector performance. In particular, HD168443, located 
within 2.5 degrees of the galactic plane, shows extremely high levels of cirrus noise (30 
mJy /aperture, approximately equal to the level predicted from the IRAS Sky Survey Atlas) 
and is not considered further. 

While we can directly examine the overall background noise in each of our 70 /im images, 
the systematic error terms are more difficult to evaluate. These errors must certainly be as 
large as those in the 24 //m data since the photospheric fitting errors should be the same for 
both wavelengths. Within our data, we find that the observed dispersion of fiuxes at 70 /im 
is consistent with larger systematic errors than at 24 fim. We assume that the systematic 
errors in the 70 /im data are 15% of the stellar fiux, about twice the dispersion in the 24 /xm 
data. Given the complex calibration required for the 70 /im Germanium photoconductors, 
this is an excellent level of performance. 

Adding all of the noise sources (photon noise, sky background, model fitting, and residual 
calibration issues) together in quadrature gives us a final noise estimate for each 70 iim target. 
In Table 2 we list these noise levels, along with the measured and the photospheric fluxes, for 
each of the planet-bearing stars. We use these noise estimates to calculate xyo, the statistical 
significance of any IR excess 

_ F70 - 

X70 = — 77 (1) 

where F70 is the measured fiux, F^, is the expected stellar fiux, and A^7o is the noise level, all 
at 70 iJ,m. The resultant histogram is shown in Fig. 4. Based on this criterion, we find that 6 
out of 25 planet-bearing stars have a 3-a or greater excess at 70 iim: HD33636 (BD+04858), 
HD50554 (BD+241451), HD52265 (HR2622), HD82943, HD117176 (70 Vir), and HD128311 
(GJ3860). Of the remaining stars, 3-a upper limits on any excess flux range from 0.1 to 1.9 
times the stellar flux, with a median of upper hmit of 0.7 F^,. The focus of this paper is the 
planet-bearing stars; the non-planet stars with excess will be discussed in a separate paper 
on the Volume Limited Sample (Bryden et al. 2005). 

Fig. 5 shows the 70 iim images for each of the six planet-bearing stars identifled as having 
signiflcant 70 /im excess. The full-width half-maximum for the telescope's point spread 
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f unction at 70 jim is 17"; there is no evidence for any emission extended beyond the point 
spread function. The background noise fluctuations are generally smooth across the field, 
with the exception of HD52265 which has obvious cirrus contamination concentrated toward 
the NW corner of the image. Within the entire ~5x2.5 arcmin MIPS field of view, most 
of the images contain point sources other than the target star. Some of these neighboring 
sources have corresponding 24 /xm emission, while others are only detected at 70 iim. 

2.3. Notes on individual sources 

HD 82943: This is our strongest detection of an excess, with a 70 /im flux 15 times the 
stellar photosphere. The lack of excess at 24 jim. requires the spectral energy distribution to 
rise steeply above the photosphere between 24 and 70 /im. As the brightest excess detection 
(113 mJy), this is the most promising candidate for follow-up detection at longer IR and sub- 
mm wavelengths. We also note that the claimed detection of ^Li in this star's atmosphere 
could be explained by the accretion of several of rock (Israehan et al. 2001, 2003). While 
our observed IR excess implies that collisions could have recently showered the star with 
some solid material, it is not clear whether the magnitude of this accretion is large enough 
to account for the enhanced lithium abundance. 

HD33636 & HD 50554: These similar stars both have considerable excess at 70 //m (~7 
times the stellar flux), but still only emit 33 and 40 mJy respectively. Although dim, these 
systems may still be detectable at sub-mm wavelength if the dust is cold enough (<, 40 
K), and even upper limits at longer wavelengths would greatly help to constrain the dust 
temperature, and hence location. 

HD128311: This star shows a deflnite detection above the photospheric level at 70 jjm, 
but a high degree of galactic cirrus and/or anomalous instrument artifacts make follow-up 
observations important for determining the strength of the excess more accurately. The 
Spitzer/IR^ spectrum for this source lies well below the low quality (FQUAL=2) detection 
at 25 /um by IRAS and shows no evidence for an excess out to 30 /xm. Additional photometric 
measurements with MIPS are currently scheduled. 

HD 117 176 (70 Vir): Unlike all of our other stars with 70 fim excess, the accuracy of the 
flux measurements for this star are limited by systematic/fltting errors, rather than by sky 
noise. Additional work to determine the source of the systematic errors will help quantify 
the overall noise level for this star. 

HD52265: Although this star has a high level of excess at 70 /^m (a factor of 4 above 
the stellar photosphere), this is our least statistically signiflcant detection of excess due to 
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a high level of galactic cirrus in the field. This noise (7.4 mJy) unfortunately may not be 
greatly reduced with additional integration time. Interestingly, this star happens to fall in 
the limited field of view of the future space-based telescope COROT (Baglin et al. 2002). As 
the only planet-bearing primary target for that mission, its visible emission will be monitored 
over an extended time period for any variability due to transiting objects. 

HDI42: This star has a fiux twice that of the stellar photosphere, but the potential 
excess is only significant at the 2.4-0" level. A relatively short integration time was used (3 
cycles) resulting in above average background noise. Additional follow-up observations of 
this star should clarify whether the excess is due to dust or is just a noise fiuctuation. 

55 (p) Cnc: With a measured flux of 18.9 ± 3.8 mJy (compared to a predicted stel- 
lar flux of 18.2 mJy), we flnd no evidence for an IR excess. This is in contrast with a 
claimed detection of 60 /im excess by ISO (Dominik et al. 1998) at the level of 170 ± 30 
mJy. The original detection has more recently been attributed to background extragalactic 
sources based on sub-mm mapping of the region (Jayawardhana et al. 2002). We conflrm 
the presence of two (presumably) extragalactic 70 //m sources located at 08:52:37-1-28:20:05 
and 08:52:35-1-28:19:18, which correspond to sub-mm sources 1 and 3 from Jayawardhana 
et al. (2002). However, the 70 /xm flux from these neighboring sources is small; a beam large 
enough to enclose them makes little appreciable change in the total integrated flux. 

p CrB & HD210277: Along with 55 Cnc, these stars were identifled as having extended 
emission via coronographic measurements at visible wavelengths (Trilling et al. 2000). We 
flnd no evidence for excess emission around either star at 24 or 70 pva.. 

3. Correlation of Excess with System Pctrameters 

Here we examine whether IR excess is correlated with stellar properties (age, metallicity, 
spectral type) and planet orbital characteristics (planet mass, semi-major axis, eccentricity). 
The planet and star properties are summarized in Tables 1 and 3. 

The planetary systems considered here all have at least one planet with mass greater 
than Jupiter's and with semi-major axis smaller than Jupiter's. Considering this sample 
in aggregate, there is a clear trend between the presence of this type of planet and the 
detection of IR excess. Within our sample we have found 6 stars with excesses out of the 
25 planet-bearing stars for a detection rate of 24 ± 10%. This is much less than the excess 
detection rates for young stars; Habing et al. (2001) observe that 60% of A-K stars younger 
than 400 Myr have excess 60 /im emission, while Rieke et al. (2004) similarly find that half of 
young A stars have excess 24 pm. emission. Analysis of IRAS and ISO data suggests that the 
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fraction of stars with excesses drops considerably with increasing age, however, particularly 
beyond 1 Gyr (Habing et al. 2001; Decin ct al. 2003). The nominal IRAS statistics suggest 
that excesses are found around just 10-15% of all main- sequence stars (Backman & Paresce 
1993), a factor of two less than our detection rate. If ages older than 1 Gyr are considered, 
the disparity in rates becomes more significant. ISO surveys found excesses around only 
10±3% of stars older than 400 Myr (Decin et al. 2000; Habing et al. 2001). This detection 
rate is similar to that found within our sample of stars not known to have planets (Bryden et 
al. 2005). Within this group, 10 ±4% (6 out of 58) of the stars exhibit excess IR emission at 
70 /im. Among stars older than 1 Gyr, the detection rate for stars with planets is more than 
double that for stars without. This tentative result would be consistent with the suggestion 
of Dominik & Decin (2003) that excesses due to Kuiper Belt material can have very long 
lifetimes. The apparent increase in frequency of dust disks around planet-bearing stars over 
stars without known planets implies that searches for dust disks might benefit from targeted 
searches of planet-bearing stars. The converse would also be true - planet searches, vis radial 
velocity or otherwise, may benefit from specifically targeting stars with IR excesses. This 
correlation is preliminary and will be examined closely as additional stars in this survey are 
observed. 



What kind of planets do stars with debris disks possess? Table 1 and Fig. 6 summarize 
the properties of the planetary systems. Two of the systems with strong excesses have large 
planets on distant, eccentric orbits (9.3 Mjup at 3.6 AU for HD33636 and 4.9 Mj^p at 2.4 AU 
for HD50554). The system with strongest excess, HD82943, has two moderate-sized planets 
with masses ~1 Mj^p at 0.73 and 1.16 AU. Overall, however, the distribution of planets 
associated with IR excess is little different than for those without an excess. For every type 
of planet with IR excess, there are counter-examples without excess: for example, HD39091 
has a large planet on a distant, eccentric orbit yet shows no dust emission. The only possible 
trend is that no systems with short-period planets (<0.25 AU) exhibit any IR excess, though 
the small-number statistics limit the strength of this conclusion. 

Planets should affect the diffuse material in their immediate vicinity. The region under 
the direct gravitational influence of a planet at with eccentricity has a half-width of 



y^4e2a2/3 + i2r|^, where rn = ap{M^/^M^f/^ is the Hill radius (e.g. Bryden et al. 2000). 

For the planet orbiting HD33636 (MpSim=9.28A'/jup, 0,^=3.56 AU, 6^=0.53) this directly 
perturbed region extends from 1 to 6 AU where we have ignored the unknown inclination 
effect which enters only weakly into the Hill radius. This zone would be cleared of any large 
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objects that might provide a local source of dust generation. Although dust could drift 
toward the planet from a distant belt of source bodies, Jupiter-mass planets are efficient 
at gravitationally scattering this dust from their neighborhood (Moro-Martin & Malhotra 
2005). Outside of the planet's region of direct influence, substantial Kuipcr Belt material 
could be present at distances ^ 10 AU and effective temperatures below 100 K, consistent 
with the cool excess (no 24 /xm emission) seen toward these stars. Additional Spitzer and 
ground-based submillimeter data are required to constrain the temperature and the physical 
distribution of this material (Sheret et al. 2004). 



3.2. Stelleir Properties 

Stellar ages cited in the literature for these stars span at least a factor of two, highlighting 
the difficulty in determining ages for mature, main sequence stars. Whenever possible, we 
adopt ages from the compilation of Wright et al. (2004) which provides a uniform tabulation 
for 1200 stars based on Ca II H&K fine strengths. While Table 3 lists the ages for the 
planet-bearing stars. Fig. 7 gives a histogram for the larger sample, both with and without 
planets. 

There is no obvious age difference between the overall planet-bearing and non-planet 
bearing samples; their median ages are 6 and 4 Gyr respectively. The one outlier among 
the planet-bearing sample is HD128311, which has no well-determined age. It is hsted as a 
possible member of the UMa moving group (King et al. 2003) based on its chromospheric 
activity and its X-ray detection in the ROSAT All Sky Survey. Both of these indicators 
suggest an age of 0.5-1 Gyr. In Fig. 7, the stars with excess are flagged with arrows, solid 
for planet-bearing stars and open for those without. With the exception of HD128311, the 
planet-bearing stars with excesses are all much older than a billion years. The planet-bearing 
excess stars are slightly younger than the planet-bearing sample as a whole (4.4 vs. 6.0 Gyr), 
but the difference is not statistically significant. Any inverse correlation of age with amount 
of dust is much weaker than that observed in stars younger than 1 Gyr (e.g. Decin et al. 
2003; Rieke et al. 2004). 

Fig. 8 presents a similar histogram for the metallicity of stars with/without planets, 
again marking IR excess stars with arrows. The distribution of metallicities for planet- 
bearing stars (dark grey) is clearly different than for those without planets (fight grey). 
The typical metallicity of stars with planets is well known to be greater than solar (Wright 
et al. 2004; Fischer & Valenti 2003) so the fact that the average metallicity of these stars 
is slightly greater than solar, [Fe/H]= 0.13 ± 0.05 with a dispersion of 0.19, is no surprise. 
Although the planet-bearing stars with excess have an average metallicity above solar, their 
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average [Fe/H]=0.05 is below the planet-bearing average. Comparing the stars with and 
without planets to their own samples, there is no evidence for higher metallicity resulting in 
a greater amount of IR emitting dust. 



4. Analysis of Detected Dust 

Unfortunately, our detections of IR excess provide only limited information about the 
properties of the dust in each system. Ideally, we would like to determine the dust's temper- 
ature, mass, grain sizes, and overall distribution. Because none of our sources are resolved 
(see Fig. 5), we have no information on whether the dust disk has any asymmetric structure 
such as a warp, a centering offset, or general dumpiness. Still, we can use the spectral en- 
ergy distribution (SED) to provide some constraints on the dust properties. The most direct 
constraint provided by the SED is on the dust temperature, which then can be converted 
into a radial distance from the central star. As seen in Fig. 3, 24 /xm measurements provide 
an upper limit on the dust temperature, while sub-mm limits set a lower limit. Because we 
usually have no information longward of 70 /im, however, only an upper limit on the temper- 
ature (or an inner limit on the dust's orbital location) can be derived. If a dust temperature 
is assumed, the observed flux can be translated into the total dust disk luminosity relative 
to its parent star. 

For disks with detections of IR excess, a minimum dust luminosity can be calculated. 
Using the blackbody formula for flux (F oc L/T^{e^^/'^'^ ~ 1))) for ^ong wavelength radiation 
(i.e. on the Rayleigh- Jeans tail of the stellar blackbody curve), the ratio of dust to stellar 
fluxes becomes 

-^dust _ -^dust hvT^ 



A;Td^^,,(eWfe^dust _ 1) 

The minimum disk luminosity as a function of 70 /xm dust flux can be obtained unambigu- 
ously by setting the emission peak at 70 //m (or, equivalently, Tdust = 52. 5K) : 

-^dust/ . . V ,^_5 /SeOOKV'^ i^70,dust /„v 

-^(mmimum) = 10 \^ — 1 — (3) 

Consistent with this formula. Fig. 9 shows constraints on Ldust/Li, and Tdust for HD82943, 
the planet-bearing star with 70 iim flux a factor of 15 above the stellar photosphere. The 
lines in this flgure are based on the 3-a limits to the observed 24 and 70 iim fluxes. Tdust 
is meant to signify the typical emitting temperature for the dust; in reality some range of 
temperatures will be found in any given system. For this star, our strongest detected excess, 
the minimum Ldust/-^* is relatively high - 1.4 x 10~^ - more than two orders of magnitude 
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above estimates for the Kuiper belt's luminosity. (A listing of minimum Ldust/-^* for the 
other planet-bearing stars is given in Table 2.) The dust temperature is well-constrained by 
the lack of 24 fim excess to be 70 K at the 3-a level (the cross-hatched region in Fig. 9) or 
^ 60 K at the 1-0" level (the solid- filled region). Dust at lower temperatures showing greater 
emission at longer wavelengths cannot be ruled out. 

Sub-mm observations are critical to constraining the dust properties and thus the overall 
mass of the material responsible the the excess in these stars. Fig. 10 shows the limits on 
the dust characteristics for a case where sub-mm data is available to help constrain the dust 
temperature on the low end. (Greaves ct al. (2004) observed HDl 17176 at 850 /im, measuring 
a flux of 2.3 ± 2.4 mJy.) Note that in order to set the most strict limits, we assume that the 
grains are small enough that their emissivity drops off as for long-wavelength radiation 
(see e.g. Draine & Lee 1984; Wyatt & Dent 2002). The sub-mm measurement does set a 
lower limit for the dust temperature, but because the 70 /im emission is not as strong as for 
the previous star, the upper limit on tempcratiire is less strict. Overall, the possible range 
of dust temperature and luminosity arc much better constrained: 1-a limits (the solid-filled 
region in Fig. 10) give Tdust ~ 20 to 100 K, Ldust/^* ~ 10"^-2 to lO'^-^. For a temperature 
of 50 K, this disk brightness corresponds to a total dust surface area of ~10^'^-10^^ cm^ or, 
equivalently, a dust mass of ^lO^^-lO^^M® in 100 /im grains. The mass contained in larger 
planetesimals is even less certain. 

For stars with no detected emission, a 3-o" upper limit on the fluxes leads to upper limits 
for Ldust/-^* of 10^^ to 10~^, assuming a dust temperature of ~50 K (again see Table 2). 
For these stars we are setting limits of ~10-100 times the level of dust in our solar system. 
Deeper integrations on these stars may result in detections of fainter excesses or improved 
upper limits, but extragalactic confusion or cirrus will set limits on the amount of detectable 
excess for these relatively distant stars. Non-planet bearing stars in the VLS are brighter 
and typically located in more cirrus-free regions, so that our Ldust/-^* limits for these objects 
tend to be more sensitive than for the planet-bearing stars reported here. 

5. Conclusions 

Spitzer/MlPS observations of 26 main-sequence FGK stars which are known to have 
planetary companions reveal 70 /im dust emission from 6 debris disks. These are the first 
examples outside the solar system of stars which have both well-confirmed planets and well- 
confirmed IR excesses. We conclude from this study that excesses at 70 fira arc at least as 
common around planet-bearing F5-K5 stars as they are around the A-F stars first detected 
by IRAS. 
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Our finding is in contrast with Greaves et al. (2004), who found no overlap between stars 
with disks and stars with planets. In order to look for dust /planet systems, they observed a 
sample of 8 planet-bearing stars in the sub-mm, with no detections. The disks identified by 
Spitzer, however, are considerably less massive than those detectable in the submillimeter. 
For example, in the case of e Eri, the ratio of 60 fim IRAS emission (1,600 mJy) to 850 fim 
emission (40 mJy) is 40:1. If this same ratio applies to the sources seen here which have 
a typical 70 //m brightness of 35 mJy, the predicted submillimeter flux density would be 
~1 mJy which is just the l-cr noise limit for the SCUBA survey. Six of the eight sub-mm 
sources are also in our target sample, allowing a direct comparison of Spitzer with SCUBA. 
In contrast with the sub-mm non-detections, five photospheres were observed with SNR>, 5, 
and one star (HD117176, 70 Vir) was found to have significant 70 /xm excess. While SCUBA 
follow-up observations of the sources detected here will be very important for constraining 
the temperature and total mass of dust, conclusions about the incidence of disks around 
planet- bearing stars require Spitzer's higher sensitivity. 

While excesses at 70 //m are relatively frequent in our survey, we find no evidence for 
excess at 24 /im around any of the planet-bearing stars. The lack of 24 /im emission in 
these systems is consistent with the previously known lack of warm material toward stars 
older than a few 100 Myr (Aumann & Probst 1991; Fajardo-Acosta et al. 2000; Laureijs 
et al. 2002). Assuming peak 24 /xm emission from 150 K dust and the observed dispersion 
{Fd/F^ < 0.07), we can set a 3-(7 upper limit on Ldust/-^* for hotter dust of 5 x 10~^ which 
is well above the 10~^ level of emission from our asteroid belt, but well below the level seen 
toward young A stars such as (3 Pic or Fomalhaut (a PsA). As described above, the lack of 
warm material is consistent with its destruction or dispersal by the gravitational effect of 
planets in the 1-5 AU region (Liou & Zook 1999). 

For stars with a single measurement of excess at 70 //m, the dust properties are not 
well constrained, but are generally consistent with Kuiper belt configurations - distances of 
several tens of AU and temperatures ~50 K. The observed dust luminosities (relative to the 
central stars'), however, are much brighter, exceeding the Kuiper belt's Ldust/-^'* by factors 
of -100. 

Within the planet-bearing sample there is little to no correlation of 70 ^m. excess with 
stellar age, metallicity, or spectral type. But there is a strong hint of a positive correlation 
between both the frequency and the magnitude of dust emission with the presence of known 
planets. As more Spitzer observations of planet-bearing and non-planet-bearing stars are 
made, as part of this survey and in other programs, these statistics will become much better 
defined. Additional radial velocity studies of stars with excesses but without presently known 
planets will test the predictive power of this correlation. 
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HD# 


name 


Mpsim (Mj„p) 


ttp (AU) 




142 


GJ4.2A 


1.36 


0.98 


0.37 


1237 


GJ 3021 


3.21 


0.49 


0.505 


13445 


Gl 86 


4 


0.11 


0.046 


17051 


L Rot, HR 810 


2.26 


0.925 


0.161 


27442 


GJ167.3 


1.28 


1.18 


0.07 


33636* 


BD+04858 


9.28 


3.56 


0.53 


39091 


GJ 9189 


10.35 


3.29 


0.62 


50554* 


BD+241451 


4.9 


2.38 


0.42 


52265* 


HR 2622 


1.13 


0.49 


0.29 


75732 


55 (p) Cnc 


0.045 


0.038 


0.174 






0.84 


0.11 


0.02 






0.21 


0.24 


0.34 






4.05 


5.9 


0.16 


82943* 


BD-112670 


0.88 


0.73 


0.54 






1.63 


1.16 


0.41 


95128 


47 Uma 


2.41 


2.1 


0.096 






0.76 


3.73 


< 0.1 


114783 


GJ 3769 


0.9 


1.2 


0.1 


117176* 


70 Vir 


7.44 


0.48 


0.4 


120136 


T Boo, GJ 527A 


3.87 


0.0462 


0.018 


128311* 


GJ 3860 


2.63 


1.06 


0.21 


134987 


23 Lib 


1.58 


0.78 


0.25 


143761 


p CrB 


1.04 


0.22 


0.04 


145675 


14 Her 


4.89 


2.85 


0.38 


160691 


p Arae, GJ 691 


0.042 


0.09 









1.7 


1.5 


0.31 






3.1 


4.17 


0.8 


168443 


GJ 4052 


7.7 


0.29 


0.529 






16.9 


2.85 


0.228 


169830 


HR 6907 


2.88 


0.81 


0.31 






4.04 


3.6 


0.33 


177830 


GJ 743.2 


1.28 


1 


0.43 


186427 


16 Cyg B 


1.69 


1.67 


0.67 


210277 


GJ 848.4 


1.28 


1.097 


0.45 


216437 


p Ind 


2.1 


2.7 


0.34 



*star with excess 70 pm emission 



Table 1: Orbital characteristics of the known planets surrounding stars observed by Spitzer. 
Also see Fig 6 for a graphical display of some of this information. 
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24 fjLm 










70 /im 








HD# 


-^MIPS 




Fmips /Fi, 


-^MIPS 


F. 


Fmips /Fi, 


SNR 


X70 




-f'dust/ 


142 


117.3 


123.3 


0.95 


26.4 ± 5.1 


13.9 


1.9 


5.6 


2.4 




<2.0xl0-^ 


1237 


82.9 


88.7 


0.94 


10.0 ± 2.9 


10.1 


1.0 


4.0 


0.0 




<9. 0x10-6 


13445 


162.9 


166.9 


0.98 


3.9 ± 6.3 


19.0 


0.2 


0.7 


-2.4 




<4.1xl0-6 


17051 


166.8 


161.7 


1.03 


20.1 ± 4.1 


18.1 


1.1 


6.4 


0.5 




<6. 7x10-6 


27442 


1147.7 


1222.7 


0.94 


126.5 ± 21.9 


139 


0.9 


19.2 


-0.6 




<7.8xl0-6 


33636* 


42.9 


44.1 


0.97 


33.1 ± 2.3 


5.0 


6.7 


15.4 


12.4 


32.4 


4.9x10-5 


39091 


139.9 


150.4 


0.93 


21.5 ± 3.6 


17.0 


1.3 


8.3 


1.3 




<9.0xl0-6 


50554* 


47.1 


51.2 


0.92 


39.5 ± 2.8 


5.8 


6.8 


14.9 


12.1 


38.8 


4.4x10-5 


52265* 


74.3 


70.7 


1.05 


31.7 ± 7.4 


8.0 


4.0 


4.3 


3.2 


27.3 


2.9x10-5 


75732 


172.8 


162.7 


1.06 


18.9 ± 4.5 


18.2 


1.0 


5.3 


0.2 




<8. 3x10-6 


82943* 


66.5 


66.5 


1.00 


113.3 ± 6.8 


7.5 


15.1 


17.0 


15.7 


122 


1.2x10-^ 


95128 


259.5 


265.9 


0.98 


29.1 ± 6.0 


30 


1.0 


7.3 


-0.1 




<4.8xl0-6 


114783 


45.2 


52.1 


0.87 


6.2 ± 2.2 


6.0 


1.0 


3.2 


0.1 




<1.4xl0-5 


117176* 


373.6 


395 


0.95 


77.4 ± 10.2 


44.8 


1.7 


10.2 


3.2 


37.6 


1.0x10-5 


120136 


338.2 


321.7 


1.05 


30.9 ± 7.9 


36.3 


0.9 


5.4 


-0.7 




<3.7xl0-6 


128311* 


67* 


73.5 


0.91 


26.5 ± 3.9 


8.4 


3.2 


7.1 


4.6 


20.8 


3.0x10-5 


134987 


78.8 


68.3 


1.15 


3.1 ± 4.5 


7.7 


0.4 


0.7 


-1.0 




<1. 1x10-5 


143761 


201.8 


192.5 


1.05 


27.8 ± 6.1 


21.7 


1.3 


5.4 


1.0 




<9. 6x10-6 


145675 


92.4 


106.4 


0.87 


9.2 ± 2.5 


12 


0.8 


5.4 


-1.1 




<4. 8x10-6 


160691 


270.6 


273.5 


0.99 


24.1 ± 8.5 


30.8 


0.8 


3.4 


-0.8 




<6. 1x10-6 


168443 


55.6 


53.3 


1.04 


9.6 ± 29.6^ 


6.0 


1.6 


0.3 


0.1 




<1.5xl0-^ 


169830 


99.0 


104.8 


0.94 


4.7 ± 5.9 


11.9 


0.4 


0.8 


-1.2 




<6. 5x10-6 


177830 


88.6 


104.3 


0.85 


4.4 ± 3.5 


11.8 


0.4 


1.5 


-2.1 




<3.2xl0-6 


186427 


93.0 


103.1 


0.90 


3.0 ± 5.6 


11.6 


0.3 


0.5 


-1.5 




<4.3xl0-6 


210277 


83.5 


91.9 


0.91 


8.0 ± 2.9 


10.4 


0.8 


3.2 


-0.8 




<5. 1x10-6 


216137 


107.5 


102.1 


1.05 


9.5 ± 3.9 


11.5 


0.8 


2.7 


-0.5 




<8. 1x10-6 



*star with excess 70 /^m emission 

Must fluxes have been color corrected by 15%, appropriate for ~50K emission 
*24 /xm flux is from IRAS 25 /im 
*high galactic cirrus contamination 



Table 2: Measured and predicted fluxes at 24 and 70 //m for the planet-bearing stars, in mJy. 
For the 70 //m data, wc also list the signal-to-noise ratio (SNR), the signiflcance level of any 
excess {xvo', see Eq 1), and L^ust/Li,, the minimum disk luminosity (see Eq 3 in §4). 
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*star with excess 70 jim. emission 

^see discussion of HD128311's age in §3.2 

Table 3: Spectral type, metallicity, distance, and age for planet-bearing stars observed by 
Spitzer. 
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F/F^, at 

Fig. 1. — Distribution of 24 //m fluxes relative to the expected photospheric values for a 
sample of 84 stars. A gaussian distribution with dispersion 0.07 is shown for comparison. 
One star (HD69830; Beichman et al. 2005) clearly stands out from the main population of 
stars with no excess emission above their stellar photospheres. In this histogram and those 
that follow, stars with planets are shaded with dark grey while stars without known planets 
are designated by light grey. 
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Fig. 2. — Distribution of 70 //m fluxes relative to the expected photospheric values for a 
sample of 84 stars. While most stars cluster around 1, where their flux is photospheric, 
several stars show a high degree of excess emission. Stars with planets are shaded as dark 
grey while stars without known planets are designated with light grey. Although planet- 
bearing stars make up less than a third of the sample, four out of the five stars with the 
highest factor of excess 70 /xm emission are known to have planets. 



-25- 



A (yLtm) 

10 



100 



1000 



-| — I — I I I I I 



-| — I I I I I I 



-| — I — I I I I I 



-i — I — I I I r 



HD82943 



Tdust=50K 




_l I I I I I 



1000 



Fig. 3. — Spectral energy distributions for HD82943 and HD 11 71 76. In addition to our 24 
and 70 /^m Spitzer data, (dark circles), we also show optical measurements and IRAS fluxes 
at 12 and 25 //m. For HDl 17176, a sub-mm (850 jim) constraint is also available (Greaves 
et al. 2004). The emission from dust at a given temperature (dashed hues) is added to the 
stellar Kurucz model (dotted line) in order to fit the observed excess emission at 70 /im. In 
each plot, two separate fits (two different dust temperatures) are considered. In the case of 
HD82943 hot dust (150 K) is ruled out by the 24 /jm observations, while for HD117176 cold 
dust (20 K) is excluded by the sub-mm upper limit (see §4 for discussion of the limits on 
dust temperature and luminosity). 
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Fig. 4. — Statistical significance of deviations from the expected stellar flux at 70 /xm. A 
gaussian distribution of errors is shown for comparison. In order to be considered a significant 
70 fjLm excess, X7o should at least exceed the 3-a level. Six of the stars with planets (dark 
grey) meet this criterion, as do six of the stars without known planets (light grey). 



Fig. 5. — 70 nm images for the six planet-bearing stars with excess emission. The contrast 
scale is the same for each image. The images arc ~2 arcmin on each side; the beam size 
fwhm is 17". North is oriented upward and east to the left within each frame. 
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Fig. 6. — Orbital characteristics for the survey stars. Minimum planet mass (MpSini) is 
plotted versus each planet's semi-major axis (a^). The lines extend from periapse (ap(l — e^)) 
to apoapse (ap(l + e^)). The planets which are associated with 70 /im excess are marked as 
triangles, while those without are circles. 
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Age (Gyr) 

Fig. 7. — Distribution of stellar ages. Stars with planets are shaded as dark grey while stars 
without known planets are designated with light grey. The ages of stars with 70 fim excess 
are flagged as arrows at the top of the plot. For planet-bearing stars the arrows are filled; 
for stars without known planets they are open. The length of the arrow is an indicator of 
the strength of 70 iim excess. 
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Fig. 8. — Distribution of stellar mctallicities. Stars with planets are shaded as dark grey 
while stars without known planets are designated with light grey. The metallicities of stars 
with 70 //.m excess are flagged as arrows at the top of the plot. For planet-bearing stars the 
arrows are filled; for stars without known planets they are open. The length of the arrow is 
an indicator of the strength of 70 //m excess. 
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Fig. 9. — Constraints on the temperature and total luminosity of the dust around HD82943, 
as provided by 24 and 70 yUm MIPS observations. The area to the upper right (hot temper- 
ature, bright luminosity) is ruled out by the 24 /im upper limit on the dust emission. The 
3-cr limits around 70 /xm detection confine the dust characteristics within the cross-hatched 
region. The sohd-fiUed region is defined by 1-a error fimits. Without any longer wavelength 
information, there is no lower limit on the dust temperature or upper limit on the dust 
luminosity. The approximate characteristics of the asteroid and Kuiper belts are shown for 
comparison. 
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Fig. 10. — Constraints on the temperature and total luminosity of the dust around 
HD117176, as provided by 24 and 70 iim MIPS data combined with sub-mm observations. 
The upper right is ruled out by the 24 /im upper limit, while the upper left is ruled out by the 
sub-mm (Greaves et al. 2004). Based on the 70 //m data, the dust temperature and luminos- 
ity are confined to the shaded region. The 3-a limits for the three wavelengths confine the 
dust characteristics within the cross-hatched region, while the solid-filled region is defined 
by 1-a error limits. Again, the approximate characteristics of the asteroid and Kuiper belts 
are shown for comparison. Note that although the formal 3-a error limits extend as low as 
the Kuiper belt's luminosity, the Kuiper belt itself would be too faint to detect. 



